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Technological applications of many-body stmctures that emerge in gated devices under minimal con¬ 
trol are largely unexplored. Here we show how emergent Wigner crystals in a semiconductor quantum 
wire can facilitate a pivotal requirement for a scalable quantum computer, namely transmitting quantum 
information encoded in spins faithfully over a distance of micrometers. The fidelity of the transmission 
is remarkably high, faster than the relevant decohering effects, independent of the details of the spatial 
charge configuration in the wire, and realizable in dilution refrigerator temperatures. The transfer can ev¬ 
idence near unitary many-body nonequilibrium dynamics hitherto unseen in a solid-state device. It could 
also be useful in spintronics as a method for pure spin current over a distance without charge movement. 


Introduction.- Spin chains can facilitate several impor¬ 
tant technological applications such as pure spin currents in 
non-itinerant systems JH], quantum state transfer lii and 
quantum gates | 0 - 0 ]- Most of the above applications are 
facilitated by a nearly unitary dynamics of the spin chain. 
Such dynamics is not only interesting for potential quan¬ 
tum technology J^, Si but also fundamentally important to 
address questions of equilibration quantum thermodynam¬ 
ics and information propagation iS]. Thus, the physical re¬ 
alization of artificial spin chains that have the potential for 
long time unitary dynamics is an important quest - they 
have been realized only very recently, and exclusively in 
atomic physics systems: in cold atom systems lUOL [Till , 
ion traps and Rydberg systems IB^ . In the realm 

of solid-state, on the other hand, nonequilibrium dynamics 
of engineered spin chains is far from unitary and primar¬ 
ily driven by equilibration/relaxation ifT^ . Although some 
bulk magnetic materials ifl^ [Tv |, NV centre chains 10, @] 
and Josephson junction arrays 111^ hold the potential for 
long chain unitary dynamics, that is still somewhat distant 
from experimental realization. As far as the semiconduc¬ 
tor realm is concerned, permanently fabricated spin chain 
structures with dangling bonds ifl^ or phosphorus dopants 
lE5l may also hold the potential, but is yet to be examined 
either theoretically or experimentally. A natural question is 
thus whether one can realize spin chains exhibiting nearly 
unitary nonequilibrium dynamics in a feasible manner in a 
two dimensional electron gas (2DEG) and thereby open the 
door to the aforementioned applications in solid state. 


Individual electrons trapped in gate defined quantum dot 
arrays have been proposed for simul ating spin chains in 
the limit of one electron in each dot 12 ih . Although the 
fabrication of larg e quantum dot arrays is an active cur¬ 
rent effort ll^l23n and the complex electronics for gate ad¬ 
dressing is also being designed lE^ . it is worth considering 
the potential of simpler gated structures. In particular, self 
assembled charge configurations, such as Wigner crystals, 
are naturally formed without demanding local control. By 


controlling the density of electrons one can vary the dis¬ 
tance between the charges and consequently engineer their 
exchange interactions. This, in comparison with quantum 
dot arrays, allows for stronger exchange couplings which 
then potentially provide more thermal stability, faster dy¬ 
namics and less sensitivity to decoherence. 

Here we show that emergent “self assembled” electronic 
spin chains arising due to Wigner Crystallization in quasi- 
ID nano-wires can be probed with two spatially separated 
accessible interfaces (two quantum dots) so that nonequi¬ 
librium dynamics and its applications can be probed. Par¬ 
ticularly we show how this setting can be used to transfer 
spin qubits between two quantum dots separated by pm 
scales, which is currently being actively considered as an 
important problem, with very few suggested solutions & 
Unn . Additionally we suggest a feasible way of observing 
this phenomenon through pA scale currents, which, in turn, 
opens up a new option in low dissipation spintronics for a 
spin current without a charge current. 

Wigner Crystalization.- By applying strong confining 
potentials on a 2DEG one may trap a few electrons in a 
quasi-ID region and effectively make a nano-wire. In such 
nano-wires, when the electron density is below a critical 
value, the Coulomb interactions between electrons over¬ 
takes their kinetic energies resulting in a quasi-ID Wigner 
crystal in which the electrons are extremely localised near 
to the classical equilibrium configurations. In quasi-lD 
nano-wires where the electrons are strongly confined in 
two directions, a Wigner crystal is predicted to emerge 
when the average electron-electron distance is greater than 
Aob, where is the Bohr radius ll^ . 

Model- We consider trapping N electrons (with N 
even) in a quasi-ID region by using surface electrodes over 
a 2DEG in GaAs. The trapping potential is modelled as 

/ (x-d/2)^ {x + d/2p 

VTR{x,y) = - 171 *^ 1 “^ + ho ie +e 

where m* is the electron effective mass, G is the strength 
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FIG. 1. (color online) (a) Transfer protocol for an even number of electrons in which hf, changes from 3meV to 50/reV instantaneously 
to induce dynamics, (b) classical ground state equilibrium positions for 10 electrons with fl = 550fieV/h. Coupling parameters (in 
the unit of ^.eV) are shown for the different 2- and 4-body exchange processes. Since the couplings are symmetric, only half of each 
coupling is given, (c) Similarly to (b), with G = AA0fj,eV/h. 


of the transverse potential and the two Gaussian potentials 
with height ho and width Wout define a nano-wire of length 
d extending between x = ±d/2. Two small quantum dots, 
for trapping single electrons, are formed at both ends of the 
wire by applying proper voltages to the gates. The quantum 
dots are modelled by the following potential 


(x + d/2~l)^ (x-d/2+lf 

VQD{x,y) = hb [e ^ -F e ^ 
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where the first two Gaussian potentials with heights hb and 
width w, centered at the distance I from the boundaries of 
the wire, create two quantum dots at both sides of the nano¬ 
wire and the potentials in the second line break the mirror 
symmetry by displacing the minima in the two quantum 
dots in opposite directions. The symmetry broken system 
will have a single equilibrium position for electrons and 
thus numerical convergence is easier to reach. Choosing 
appropriate values of the parameters xo, yo and a is dis¬ 
cussed in the supplemental material. 

We consider = 10 electrons trapped over a distance 
d = 1.25/im, such that the quantum dots confine one elec¬ 
tron each and the remaining 8 electrons are placed in the 
wire as shown in Fig.dla). The whole potential is 


the two endmost electrons can act as sender and receiver 
when sending quantum information through the chain. Un¬ 
like the middle electrons, we assume full control over these 
two end electrons for initialization and measurement. 

Exchange couplings - At low temperatures the dynam¬ 
ics of the Wigner crystal is approximately gov erned by the 
multi-spin-exchange (MSB) Hamiltonian 1^ . To calcu¬ 
late the exchange couplings in this Hamiltonian we exploit 
the semi-classical path integral instanton method 13014^ . 
This approach is applicable for the regimes where the quan¬ 
tum effects are small, like the case for Wigner crystals 
where the electrons are well separated (see the sup¬ 
plemental material for more details). While our methodol¬ 
ogy assumes a closed system it is potentially amenable to 
extension to directly incorporate an environment following 
the methodology of Ref. 

Using the instanton method, it was found that only pro¬ 
cesses involving up to 4*-nearest neighbour pairwise ex¬ 
change and up to 4-body exchange were significant. With 
this restriction, the full spin Hamiltonian can be written 
(see supplemental material for the details) 
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The barrier hb is varied from 3meV to 50peV to decouple 
or couple the quantum dots to the wire, respectively. Thus, 


where cr„ = (cr^, cr^, cr*) is the Pauli vector acting on site 
n and E := (^CTj ■ cri-^i^cri ■ CTm) “F {(Tj ■ cr^)(crj. ■ 
{(jj-(ji){cTk-<ym)- Exchange couplings for U = 550p,eV/h 
and U = AAOpeV/h are given in units of peV in the ta¬ 
bles next to each charge configurations of Figs. [T](b) and 
(c). The general features are low couplings at the boundary 
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(due to the barrier between the dots and the chain) and a U- 
shaped coupling along the chain. The next-nearest neigh¬ 
bour couplings are always ferromagnetic (i.e. ,n+2 < 0) 

for the linear configuration while for the zig-zag geometry 
they show a more complex pattern varying from negative 
to positive values. The U-shaped feature of the couplings 
is because the off-center electrons are pushed towards the 
boundaries due to an unbalanced Coulomb repulsion from 
the majority of electrons on the opposite side. This makes 
the effective distance between the electrons shorter in the 
boundaries and thus results in stronger couplings. 

Quantum communication.- We assume that the quantum 
dots are initially decoupled from the wire (i.e. hf, is large). 
Furthermore, we consider zero temperature so that the elec¬ 
trons are confined to their lowest vibrational mode and are 
prepared in their spin ground state IGS"). The electron in 
one quantum dot is prepared in an arbitrary quantum spin 
state lip) = cos(0/2) |t) + sin 0/2 |j,), which is sup¬ 
posed to be transferred to the opposite dot, which confines 
an electron in unpolarized mixed state 1/2. The initial state 
of the system is thus/9(0) = |'!/’)(V'|C>|GS')(GS'|(E>//2. The 
barriers are then simultaneously lowered (i.e. hb is low¬ 
ered) in order to start a unitary dynamics under the action 
of the new Hamiltonian. One can compute the density ma¬ 
trix of the last electron pN{t) by tracing out the others. To 
quantify the quality of transfer one can compute the fidelity 
F{t) = {'il!\pN{t)\ip) which is independent of 6 and cp due 
to the SU(2) symmetry of the Hamiltonian. 

In Figs. 12 a) and (b) we plot F{t) as a function of time 
for the two charge configurations with Q = 550/ieU/fi 
(linear chain) and 0 = AAOpeV/h (zig-zag chain). It is 
clear from these figures that the fidelity peaks at a time 
t = tm and the takes its maximum value Fmax = F{tm)- 
Although the linear chain gives a faster dynamics with 
tm — 4ns (due to fairly larger couplings) in compari¬ 
son with the zig-zag charge configuration with tm — 8ns, 
the maximum fidelity Fmax is remarkably high for both 
configurations, certainly larger than a uniform chain IISTII . 
These results illustrate the key point that the details of 
the charge configurations are not important for the qual¬ 
ity of spin transport. To see the scalability, we plot Fmax 
as a function of N in Fig. |2c), keeping the density of 
electrons fixed and using only nearest-neighbour couplings 
since these are by far the largest and computing higher or¬ 
der interactions becomes computationally prohibitive for 
> 10 electrons. As the figures shows the fidelity Fmax 
remains very high even for chains up to = 20 electrons. 

Explanation.- To understand the remarkably high fideli¬ 
ties achieved through the time evolution of the Wigner 
crystal, one has to consider the nearest-neighbor couplings, 
which dominate the Hamiltonian FI. As shown in the ta¬ 
bles in Fig. [T] (b) and (c), the couplings J2,3 (which are 
identical to Jn- 2 ,n-i) are almost four times stronger than 
the couplings Ji,2 (which are identical to Jn-i,n) in both 
charge configurations. This implies that the ground state 
and first excited state of the Hamiltonian H show de- 





FIG. 2. (color online) Average fidelity Fav as a function of 
time, when hb suddenly changes from 3meV to SOpeV, using 
both real and effective Hamiltonians for: (a) the linear chain 
(U = 550peV/hy, and (b) the zig-zag chain (U = iAOpeV/h). 
(c) Fmax versus N. (d) The variation in Fmax when a random 
magnetic field of variance Bnuc is included in the Hamiltonian. 


localized strong correlations between the boundary elec¬ 
trons. By computing the reduced density matrix of the 
two ending spins from the ground state of H one 
can see that {'ip~\pi^]^\'ip~) > 0.8, where \'ip~) is the sin¬ 
glet state, for both charge configurations. Similarly, by 
computing puv from the first excited state of H one gets 
{'ip~^\pi^N\tp~^) > 0.9, where |'i/+) is the triplet state. The 
delocalized eigenvectors of the Hamiltonian create an ef¬ 
fective RKKY-like interaction between the two boundary 
electrons in the dots namely H^f f = Jef fcri -ctn 

where J^jj = AE /4 in which AE stands for the energy 
gap of the Hamiltonian H. In Figs.[2a)-(b) the time evolu¬ 
tion of average fidelities using both / and H are plot¬ 
ted which shows that the effective Hamiltonian can qual¬ 
itatively explain the real dynamics of the system. In fact, 
the effective model becomes more precise by decreasing 
the boundary couplings J 12 and 
Imperfections.- So far, we have assumed that the sys¬ 
tem operates at zero temperature and thus the electrons in 
the wire are initialized in their ground state. In order to 
guarantee that the proposed protocol remains valid at finite 
temperature T, one has to satisfy < AE, where ks 
is the Boltzmann constant. For the given set of couplings 
the energy gaps are AE ~ 0.77 peV for the linear and 
AE ~ 0.45 peV for the zig-zag configurations giving the 
range of temperature as ksT ^5 — 10 mK which can be 
achieved in current dilution refrigerators 
In GaAs hetero-structures the electron spins interact with 
the nuclear spins of the host material. Due to the very 
slow dynamics of nuclei spins in comparison to the time 
scales of our protocol one can describe their average effect 
on electron spin n as an effective random magnetic field 
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FIG. 3. (color online) Schematic diagram of a possible physical 
implementation of the mechanism in GaAs. Polarised electrons 
(green circles) are pumped into the left quantum dot, whilst un¬ 
polarised electrons (red circles) are pumped into the right dot. In¬ 
teraction with the confined electrons in the wire (orange circles) 
swaps the electron spins and thus exchanges the polarization of 
the currents. A spin polarisation measurement on the output cur¬ 
rent can detect the transfer of polarised electrons. 


B„. While the direction of this field is fully random its 
amplitude has a Gaussian distribution Iddll 

P(B) = exp[-B • B/2i3Lc]/(2^^Lc)'/^ (2) 

in which is the variance of the distribution. The total 


Hamiltonian thus changes as 


H 




H + "S Bn-tTk- 


(3) 


In Fig. |2ld) we plot the maximum fidelity Fmax versus 
Bnuc which shows the destructive effect of the hyperfine 
interaction. The linear configuration performs better for 
larger values of Bnuc since the faster dynamics reduces the 
time exposed to nuclear spins. A realistic value for Bnuc is 
2 —6mT (i.e. ^ 0.07 —0.23peV) at which the fidelity 
of Fmax — 0.8 — 0.95 is attainable. Using spin-orbit cou¬ 
pling one may effectively suppress Bnuc to 35 neV J^. 
At experimental temperatures ('^ 50mK) thermal 


phonons are absent in the material 04411 . Moreover, the bulk 
phonon wavelengths pm) exceed electron wavelengths 
20 nm) so much that the electrons do not couple to them 
either. However, sudden quench in hb may cause the elec¬ 
trons to jiggle around their positions causing fluctuations in 
exchange interactions. As shown in the supplemental ma¬ 
terial, these high frequency vibrations (/ ~ lOOGHz) are 
integrated out during the transmission time tm — 5ns. 

Experimental Realization - The quantum dot-quantum 
wire-quantum dot system can be realized in a III-V ma¬ 
terial such as GaAs/AlGaAs heterostructure lld^ll . using 
gates and quantum point contacts technique, as shown in 
Fig. El We assume that the quantum wire has achieved 
a state of Wigner crystallisation Spin initializa¬ 

tion in the left dot and read-out in the right dot after the 


transfer time, can be performed using demonstrated tech¬ 
niques such as singlet-triplet charge measurements 

by appending double-dots to both ends of the wire. How¬ 
ever, here we show how transport measurements can also 
probe our mechanism by effectively detecting the transfer 
of a spin current. The prototype of this measurement con¬ 
sists of a left electron pump injecting polarised 

electrons (say, from a standard source 0^ ImIi I whereas 
the right pump injecting unpolarised electrons in the re¬ 
spective quantum dots as shown in Fig. El As the elec¬ 
trons are injected into the dots they interact with the elec¬ 
trons in the quantum wire as the barrier height is reduced 
and their quantum state is swapped. Thus upon exiting 
the dots, the initially polarized electrons will be unpolar¬ 
ized and vice-versa. As the interaction time for spin swap 
is about tm ^ 4ns (tm ^ 8ns) for the linear (zig-zag) 
configuration the electrons must stay in the dot for such 
time. This implies that the pumps have to inject current of 
20pA — 40pA depending on the configuration of electrons 
in the quantum wire. A spin polarization of the current car¬ 
ried off from the right dot by the right pump could then 
be detected by the halving of the conductance through a 
quantum point contact in series with the right dot. 

Discussion and Conclusions.— We have shown that en¬ 
tirely within the remit of gate defined structures, one can 
use a Wigner crystal in a quantum wire for transferring 
spin qubits between quantum dots separated by microme¬ 
ters reaching very high fidelities. Compared to varied pro¬ 
posals for mediating between separated quantum dot spin 
qubits l&El, the proposed mechanism can have various 
advantages. For example, it could have higher speed mak¬ 
ing it more resilient to decoherence, simpler fabrication 
compared to hybrid systems, and a greater spatial extension 
over few electron quantum dot mediators. Furthermore, 
thanks to the self-assembled nature of the Wigner crystal, 
our scheme does not require complex electronics to artifi¬ 
cially create a regular structure of exchange coupled spins. 
Our analysis shows the proposed protocol can operate in 
current dilution refrigerators, and remarkably the details 
of the charge configuration does not influence the perfor¬ 
mance of the system, adding additional robustness to the 
device fabrication. It also opens an avenue in spintronics 
by facilitating the spatial transport of a spin current without 
a charge current, and provides a smoking gun for the near 
unitary nonequilibrium dynamics in many-body solid-state 
systems. Alternative realizations of Wigner crystals, in liq¬ 
uid helium Js^l, ion traps El and carbon nanotubes ll^ 
can also be used for implementing our mechanism. 
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Supplemental Materials 

1. TRAPPING POTENTIAL 


The trapping potential used in this study is 

Wr*) + l"«D(rO + E4„|/,_r I • 

j<k 

where VrRirk), Vqd{yi~) are defined in the main text. An illustration of this potential is given in Fig.lS21 

The external harrier heights in VTR{Yk) were set to he /iq = 3meV and the length of the system set to d = 1.25pm 
as discussed in the main text. We choose a value of po such that the the degeneracy between the two possible ground 
states ipi and V ’2 is lifted sufficiently with respect to the temperature, where tpi is the mirror inversion of 'ip 2 about the 
x-axis. The probability of excitation from tpi to V ’2 at temperature T in the presence of the symmetry breaking potential 
— yo)^ is roughly e ~ exp (^—m*A'^{2yo)'^/2kBT). We assume that the system is realised in GaAs, at a 
temperature of T ~ 5mK, with A = 500peV//i. Inserting these parameters e ^ exp( —Vo = lOnmis used to 
give a sufficiently low error and since we would expect this level of precision to be possible in experiments. 

The internal barrier height hi, was chosen so that the coupling between the dot electrons and the chain was significant, 
but the change in the positions of the electrons was still small. A change in hb from 3meV to 50peV was found to change 
the classical equilibrium electron positions by around 0.1% of average electron separation, whilst still giving a reasonable 
coupling. The variation in dot-chain coupling for differing barrier heights can be seen in Fig. [ST] 


N 


''(R) = E 


fc=l 



FIG. SI. Variation of the dot-chain coupling J 12 as the barrier height hb is changed. The solid red line indicates hb = 50peV. 


The remaining parameters are chosen based on the electron configuration without the dot potentials present (i.e. when 
hb = A = 0). For each particular situation, the minimum energy configuration R = {xi,yi, X 2 ,y 2 , was 

found without the dot potentials, and then the parameters were chosen such that the internal barriers would sit roughly 
halfway between the first two electrons, and such that the internal barriers have as little effect as possible on the separation 
between the dot electron and the chain. Thus I = dj2 — {x\ -|- 0:2 )/2, and empirically w = I /8, Wout = I /8 was found to 
have little effect on the electron equilibrium positions. We set cr = 1/2, so that the effect of the two quadratic potentials 
for the dots is confined to within the barriers. We set A = 2meV /h so that the confinement of the dots is much stronger 
than the confinement in the wire, and to further reduce the susceptibility to thermal noise. The final exchange couplings 
that we have calculated for this model are well within the experimentally achievable regimes which justify the choice 
of above parameters. 


2. PATH-INTEGRAL INSTANTON METHOD 

Before discussing the path-i nteg ral instanton method, it will be useful to introduce a characteristic length scale tq, to 
ease the notation. Following yj], we define this as the length scale ro = (e^/47rem*Og)^/^ at which the Coulomb 
repulsion between electrons is comparable to the potential energy of the parabolic confinement in the y-direction for a 
particular strength of the y-potential chosen as Oq = 0.5meV//i based on the energy spacing in gated quantum 
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FIG. S2. An illustration of the potential V (R) described in the text. 


dots 114ih . It is also useful to define the dimensionless parameter m = r^jaB, where as is the Bohr radius, which gives 
rn ~ 10 for the chosen value of flo- 

For a system of N electrons constrained to move in 2D, we define collecfive spafial and spin coordinafes R = 
(ri,r 2 , ...,rAr) = (xi, yi, a; 2 , y 2 ---, 2/w), cr = (cri, cr 2 ,..., cr^), and we use fo denofe the element of 

R, such that R 2 j-i = Xj, i? 2 j = Uj- In the Wigner regime these electrons are localised near to N lattice sites 


R = ...,XN,yN)- 

We assume that the conditions are such that the low-energy Hamiltonian has the form in (IS7ll . Then using the path 
integral instanton approximation, the exchange coupling corresponding to a permutation P of the electron labels is ap¬ 
proximately given hy 


* = (SI) 

where Rci is the path of least action connecting the configurations R and PR, Ai = 2 for two body exchange and 1 
otherwise, and r 7 [R] is the dimensionless Euclidean action defined as 


rim ■■= 


rP 


dr 


E 


2 ^ hy/r^ 


E(R) 


P[Rc;] is fhe resulf of quadrafic flucfuafions abouf the classical path, and is given as 


Fm,,] := 


det(-a2l + ^H(0)(r)) 
^ det'(-52l + ^H(r)) 


(S2) 


(S3) 


where H is the time-dependent Hessian matrix, defined in ferms of fhe pofenfial V (R) as 




a2E(R) 

dRi{t)dRj{t) 


R=Rci 


(S4) 


In order fo numerically compute fhe exchange couplings Jp v/e firsl find fhe “classical” configuration fhaf minimises 
fhe potential energy of fhe elecfrons. The elections are initialised evenly spaced along fhe x-axis wifh a slighf zig-zag 
perturbation in fhe y-direclion. The configuration fhaf minimises fhe potential energy is labelled as R. Typical minimum 
energy configurations in a nano-wire of lengfh d = 12ro for = 10 elections are shown for fwo differenl tiansverse 
potential 0. = 550peV and Q = ddO^eV in Fig. 1 (b) and (c) in fhe main manuscripf. For large values of Q (i.e. stiong 
confinemenfs), the electrons are constrained to lie along the x-axis, whereas for smaller Q (i.e. weak confinements) a 
zig-zag pattern forms. 
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Once the minimum energy classical configuration R has been found, the classical path of least action for each permuta¬ 
tion can he numerically found hy discredsing the exchange path into M equal dme steps At = /M, so that the action 

integral becomes a Riemann sum: 


p[R] 




1 


R(R) 


1 (rj(Tm) - rj(Tm-i))^ 

2^ 2^ l^T 2 

m=l j 


Ar 


F(R™) 


(S5) 


where = mAr and R^ = (ri(rm), r 2 (Tm), ?' 2 w('rm))- This can be minimised using a steepest descent algorithm, 

as this is equivalent to the expression for the energy of a chain of beads linked together by springs. 

For our calculations, we used = 30 and M = 70, based on the work in ll^ that shows reasonable convergence with 
these parameters for 27 electrons in a Wigner crystal. The accuracy of discredsing dme in this way has also been studied 
in by comparing numerical results to an analytically solvable system. Here it is found that for M = 64, T-r = 30 the 
error is on the order of 1 %. 

Once the classical path Rc« has been found, the prefactor R[Rc;] is calculated by diagonalising the numerator and 
denominator in eqn. (IS3I) so that the determinant can be found. This matrix is calculated along the classical path, and in 
discredsed form becomes a 2N{M -f 1) x 2N{M + 1) matrix. For example, the discredsed form of the numerator in eqn. 
(IS3b (excluding the square root) is a matrix A with elements given by 




m2 


1 


rp2 ( T ^k,l — l T ^ij “f {^k)^kl 


(S6) 


3. SPIN OPERATORS FROM PERMUTATIONS 

The multi-spin exchange (MSB) Hamiltonian is given by 

where P'^ denotes a permutation operator permudng only the spins of the particles according to the permutadon P of the 
electron labels, mp is the number of 2-body swaps that P can be decomposed into. Jp is the exchange energy, which 
depends on the energy spliding between the eigenstates of a permutation P'^. 



FIG. S3. An illustration of the types of permutation processes that we consider in this work. Jmn indicates a pairwise permutation of 
electrons m and n. indicates 3- 4- 5- and 6-body ring permutations starting clockwise from electron n. 


In this appendix we detail the conversion of the MSB Hamiltonian in (IS7b into a Hamiltonian with Pauli exchange op¬ 
erators. Bor the range of parameters considered in this study, permutations of more than 5 electrons and more than nearest 
neighbour pairwise permutations were found to be 1000 times smaller than the nearest-neighbour exchange couplings over 
the regime considered here, and so are left out. 
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With this restriction, the Hamiltonian used to simulate the evolution is 


H = P'^ - 

/ ^ ’^mn mn / ^ 

mr^n l,m,neA 

l<m<n 


J, 


{R3)- 


Iran 


E 

/c,Z,m,n 6 D 

k<l<n<7n 


J, 


(fl4) ^ 


klmn 


E 


j(R5) pa 
^ j ^ jklmn'’ 


(S 8 ) 


j,k,l,m,nEC) 

j<k<n<l<m 


where e.g. means a cycle of length 4 which permutes the spins according to ^ a/, < 7 / ^ ^ ^ 

CTfe. We use the notation m ^ nto mean m and n are up to 2"‘^-nearest neighbours, l,m,n £ A means that I, m, n are on 
the vertices of a triangle, k,l,m,n £ □ means that k, I, m, n are on the vertices of a parallelogram and j, k,l,m,n £ O 
means that j, k, I, m, n are on the vertices of trapezium (see Fig. lSSb . 


The spin permutation operators can he written in terms of pairwise spin exchange interactions, using the following 
identities 115 811 : 


per 

^12 




Pl23 + P 32 I 

-^1234 + E 432 I = 

-^12345 + -^^321 


'12 ), 

pr 


\{Pl%P2^4 


'12 “ 1 “ P 23 + Pl3 ~ 

per p(T I po- per po- per , per 
'^12-^34 “T '^13-^24 ^13 

— p<^ 1 

~ 2^'^12^34 T -'^14^23 


p<^ — 1 
-*^24 ■'■I 


per 

12-^45 


- P?2Pl3 - P 14 P 23 ) + \{P{ 

^ ! pe’' pee _i_ per per per per \ ^ 1 p' 

2 V^25-f^34 + -^24-^35 “ ^25-^341 “ 2 

1 


A"3^2"4) + ^(^"5^2' 


jer per , per per per per \ 

15-'^23 ^ -^12-^35 -^13-^25^ 

5 er per , per per per per \ 

15-^34 “T -f^l3-^45 “ '^14-^35^ 


2iA"2+'Pr5 

^(A "3 + P33 + P23 + P24 + P14) - 2 ■ 


- P23 


-PS 4 


-P43) 


(S9) 

(510) 

(511) 


(S12) 


By applying the above identities to Hamiltonian ( IS 8 I 1 . this gives the overall Hamiltonian, ignoring the factors of identity 
since they do not affect the dynamics: 


H 

(y / ^ ^mn 


mn^rn ^ n 


4 E 


(fl3) 


■ ^r\ 


l<m<n 


— jj. ^[(Tk-Cri+<Tk-Crm +cr f; ■ (Tn + (Ti ■ (Tm + cri ■ (Tn + (Tm ■ CTn + P 


klmnl 


fc,/,m,n6n 

k<l<n<m 


Q ^ ^ Jj jklm H“ ^jkln H“ jrLkra~^ jnlm H“ knlm “1“ ' ^k ^k ' 


j,k,l,m,nE^ 

j<k<n<l<m 

+ ^771 ■ + CTj ■ cr^ + (Tj ■ CT; + (Tk ■ (Tm P <^k ' n + j ■ ^m] 


(S13) 


where TjUm ■= (cTj ■ crfc)(cr; ■ cr^) + {cTj ■ CTm){crk ■ cri) - {cTj ■ cri){crk ■ (Tm)- Note the minus sign in front of the 
j(^J3) jiR5) terms since they involve permutations of an odd number of particles. Overall then, this Hamiltonian can 
be written 


4 N-r 

'f-EE Jn,n-\-r^n ' ^n-\-r H“ ^ ^ Jjklrn^jklm 

j<k<l<m 


r=l n=l 


(S14) 
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where 


*^1,71+1 

Jn,n+2 

Jn^n+'i 

Jn,n+l,n+2,n+?, 

Jn,n+l,n+2,n+A 

'^ 7 , 77 + 1,77 + 3 , 77+4 
>^ 7 , 77 + 2,77 + 3 , 77+4 



(515) 

(516) 

(517) 

(518) 

(519) 

(520) 

(521) 

(522) 


4. TWO-POINT STEEPEST DESCENT ALGORITHM 

The two-point steepest descent algorithm proceeds as follows, with the positions of all of the particles at the iteration 
labelled as R„ @]: 

• Make an initial guess Rq. 

• For each step, R„+i = R„ — a„G„, where G„ = Vrp[R„], ao is a small step size (we use oq = 0.001) and 
Q!„ = G„ ■ (R.„ - R„_i)/||G„|p for n > 0 111. 

• Stop if II Gti II < Ec (convergence to a minimum), or if || R 77 1| < Sp (convergence to a point). 

Here || ■ || indicates the Euclidean norm, and the tolerance used is Ec = £p = 1 x 10“® (although in this study the algorithm 
always converged to a minimum rather than a point). A plot showing the convergence of p is shown in Fig. lS4l for all of 
the exchange processes starting at electron 1 for = 550peV. The plot shows the change in action relative to the final 
value min(p), given as (p — min(r 7 ))/min( 77 ). Clearly the 3, 4 and 5 body processes take much longer to converge than 
2-body processes, which is intuitive given the electrons tend to move less for the 2-body processes. The small spikes in 
the plot can be explained by the algorithm occasionally overshooting the local minimum. 


5. PHONON ELUCTUATIONS 


As the barrier is dropped to couple the quantum dot to the nanowire, there may be phonon oscillations induced by the 
slight change in equilibrium positions. One effect of this would be to cause the exchange couplings to oscillate. To estimate 
the effects of this on information transfer, we calculated the average change in the nearest-neighbour exchange coupling 
for different barrier heights. We found that for a change in internal barrier height from hb = 3meV to hb = SOpeV used 
in the study, the exchange couplings changed by less than 5%. Additionally, we calculated the average classical vibration 

energy = J^ k k electron where 


2 

^x,k 


2 d'^V 


m dx^. 


CO. 




2 d^V 
m dyl 


R=R 


R=R 


(S23) 
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FIG. S4. An example of convergence of the dimensionless action 77 during the 2-point steepest descent algorithm. The trapping potential 
is G = 550/reV, and all processes start from electron 1. 


These vibration frequencies are typically on the order of l-50o = 0.5 — 2.5meV for the regimes covered in this paper, 
roughly 100 times larger than the exchange coupling. We can approximately model the effects of these phononic vibrations 
as adding small random oscillations in the exchange couplings according to 

J —> J[1 + 6 sin(d;f -|- cj))] (S24) 

where cj) is a random phase and 5 is a random amplitude chosen from a normal distribution with mean 0 and standard 
deviation 0.05. Our numerical simulations show that such fast fluctuation cancels out in time evolution and its impact on 
the fidelity is negligible. This makes the proposed protocol resistive against phonon excitations created through changing 
the barrier height hf,. 











